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Abstract. We study Hopf-Andronov bifurcations in a class of random dif- 
ferential equations (RDEs) with bounded noise. We observe that when an 
ordinary differential equation that undergoes a Hopf bifurcation is subjected 
to bounded noise then the bifurcation that occurs involves a discontinuous 
change in the Minimal Forward Invariant set. 

Dedicated to John Guckenheimer; for his many contributions to our field. 

1. Introduction 

We will consider Hopf-Andronov bifurcations in a class of random differential 
equations (RDEs) 

(1) x = fx(x)+e£ t 

as the parameter A G R is varied. Here x will belong to the plane and £ t will be 
a realization of some noise. We are interested in bounded noise in which ^ takes 
values in a closed disk A Cl 2 . More specifically we will consider Hopf-Andronov 
bifurcations with radially symmetric noise where £t takes values in the unit disk. 
The RDEs without noise, for e = 0, unfold a Hopf-Andronov bifurcation as A varies. 
We denote by U the collection of all possible realizations of the noise. We assume 
that / and U are sufficiently well-behaved that the equations (1) uniquely defines 
a flow for all realizations £ of the noise. 

For the general framework of random dynamical systems we refer the reader to 
L. Arnold's book [1] (see also [5, 11]). A distinctive feature of dynamical systems 
with bounded noise is that they may admit more than one stationary measure. 
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For discrete time Markov processes these measures were studied by Doob [7], who 
showed that their supports are precisely the minimal forward invariant (MFI) sets. 
In [9] the authors adapted Doob's proof to the case of the continuous time (not 
necessarily Markov) processes generated by RDEs on compact manifolds. It was 
observed ([9, 14]) that these stationary measures can experience dramatic changes, 
such as a change in the number of stationary measures or a discontinuous change 
in one of the supports of densities. We refer to such changes as hard bifurcations. 
Given the one-to-one correspondence between stationary measures and the MFI sets 
on which they are supported, in order to study hard bifurcations, it is sufficient to 
study the bifurcations of MFI sets themselves. 

We adopt from [10] the following assumptions on (1) and its flow: 

HI. The set A is a closed disk with smooth boundary. For each x the map 
A — >• T X X given by £ H> f(x,£) is a diffeomorphism with a strictly convex image 
/(*, A). 

H2. There exist ro > and t\ > such that 

In [10], under these noise conditions the authors provided a complete classification 
of bounded noise co-dimension one hard bifurcations in phase space dimensions 1 
and 2. We call a set F C X forward invariant if 

(2) *i(F,U) C F 

for all t G M + . Denote by T the collection of forward invariant sets. There is a 
partial ordering on T by inclusion, i.e. E ■< F if E C F. We call E C T a minimal 
forward invariant (MFI) set if it is minimal with respect to the partial ordering 
^. In this context, MFI sets were shown to exist in [9]. It follows easily from the 
definitions that an MFI set for (1) is open and connected and that the closures of 
distinct MFI sets are disjoint. 

A natural assumption is further that we are given a O 1 - invariant, ergodic proba- 
bility measure P on U with the following hypothesis: 

H3. There exist ti > so that $^(av)*P is absolutely continuous w.r.t. a Rie- 
mannian measure m on X for all t>ti and all x G X. 

Under this assumption, the closures of MFIs are supports of stationary densities 
[9]. We note that the concept of MFI set is the same as the concept of positively 
invariant maximal control sets defined in the context of control theory [6] . We now 
define bifurcation of MFI sets. 

Definition 1.1. A bifurcation of MFI sets is said to occur in a parameterized family 
of random differential equations if either: 

Bl The number of MFI sets changes. 

B2 A MFI set changes discontinuously with respect to the Hausdorff metric. 

In a supercritical Hopf bifurcation taking place in (1) for e = 0, a stable limit 
cycle appears in the bifurcation at A = 0. For a fixed negative value of A, the 
differential equations without noise possess a stable equilibrium and the RDE with 
small noise has an MFI which is a disk around the equilibrium. Likewise, at a 
fixed positive value of A for which (1) without noise possesses a stable limit cycle, 
small noise will give an annulus as MFI. A bifurcation of stationary measures takes 
place when varying A. We will prove the following bifurcation scenario for small 
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£ > 0: the RDE (1) undergoes a hard bifurcation of type B2 in which a globally 
attracting MFI set changes discontinuously, by suddenly developing a "hole" . This 
hard bifurcation takes place at a delayed parameter value A = 0(e 2 ^ 3 ). We include 
a brief discussion of attracting random cycles. Finally we demonstrate the results 
with numerical experiments on a radially symmetric Hopf bifurcation with bounded 
noise. 

For recent studies of Hopf bifurcations in stochastic differential equations (SDEs) 
we refer to [4, 3, 2, 13]. In such systems there is a unique stationary measure, with 
support equal to the entire state space. Bifurcations of supports of stationary 
measures, as arising in RDEs with bounded noise, do not occur in the context of 
SDEs. 

2. Random perturbations of a planar Hopf-Andronov bifurcation 
Consider a smooth family of planar random differential equations 

(3) (x,y) = f x (x,y) + e(u,v) 

where A G M is a parameter and u, v are noise terms from A = {u 2 + v 2 < 1}, repre- 
senting radially symmetric noise. Hypothesis HI is therefore fulfilled, we consider 
noise such that also H2 is satisfied. 

We assume that without the noise terms, i.e. for e = 0, the family of differ- 
ential equations unfolds a supercritical Hopf bifurcation at A = [12]. The Hopf 
bifurcation in planar RDEs with small bounded noise is described in the following 
result. 

Theorem 2.1. Consider a family of RDEs (3) depending on one parameter X, that 
unfolds, when e = 0, a supercritical Hopf bifurcation at A = 0. For small e > 
and A near 0, there is a unique MFI set E\. There is a single hard bifurcation 
at Abif = 0(e 2 / 3 ) as e I 0. At A = Abif the MFI set E\ changes from a set 
diffeomorphic to a disk for A < Abif to a set diffeomorphic to an annulus for A > Abif . 
At Abif the inner radius of this annulus is r* = 0(e 1 ^ 3 ). 

Proof We first bring the system without the noise terms, 

(4) (x,y) =f\(x,y), 

into normal form. Note that a smooth coordinate transformation h\ and a multi- 
plication by a smooth positive function a\ (which is equivalent to a time reparam- 
eterization) , change (4) into 

(5) (x,y) = ot\(x,y)(h\)*f\(x,y) = a x (x,y)Dh(h~ 1 \x,y))f x {h~ l (x,y)). 

A smooth normal form transformation consists of a smooth coordinate transfor- 
mation h\, a time reparameterization, and an additional reparametrization of the 
parameter A. There exists a smooth normal form transformation that changes (4) 
into 

<«> (j)-(i a 1 ) + * 4 >- 

see [12] (here r = ^ x 2 + y 2 ). Applying this normal form transformation to (3) 
yields 
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where A(x, y. A) = A) (A) + 0(r) is a two by two matrix that depends smoothly on 
its arguments. Note that expression (5) implies that the noise terms (u, v) are the 
same as in (3); they only get multiplied by a matrix A(x,y, A). Changing to polar 
coordinates we find 

r = Ar-r 3 + 0(r 4 )+ £ ^ / u A ( C ° s9 



v J V sin 



r \ \ v J \ cost 
We perform a rescaling in e by putting r = e x l 3 f and A = e 2 ^ 3 X. This brings 

r 2/3 



r = e 

(8) 



£ 2/3 



£" 2 / 3 + -_ 



u \ I -sm( 
v J ' V cos^ 



were A = A (A) + £ 1 / s O(r). Multiplying by a factor e 2 / 3 and taking the limit 
e I gives 

r = Af — f 3 + (A(u, v), (cos^, sin^)) 

for the radial component. 

Noting that (u, v) lies in a unit disk, A(u, v) lies in an ellipse. If we let a and b 
be the major and minor axes of this ellipse, then by a rotation of coordinates we 
may transform the last equation into 

(9) f = Af — f 3 + au cos + bv sin 0. 

In [10] the authors showed that the boundary of an MFI set consists of solutions 
of the extremal vector fields defined by the bounded noise differential equations. 
Observe that for e = 0, (8) reads 6 = 1 and f = and its right hand varies continu- 
ously with e. Hence for e small, the r components of the extremal vector fields are 
approximately the extremal values of the f component. A simple calculation shows 
that the extremal values of the last two terms in (9) are given by: 



± Va 2 cos 2 + b 2 sin 2 6. 

Thus the boundaries of all solutions of (9) are given by the solution r± of the 
equations: 



' sm 2 I 



(10) r± = \r± -r 3 ± ± V a 2 cos 2 + b 2 j 

Now note that is of order e~ 2 ^ 3 . This allows us to average equations (10) to obtain 
the averaged equations: 

(11) f± = Af± - f 3 ± ± / Va 2 cos 2 + b 2 sin 2 OdO. 

Jo 

Note that the integral may be transformed into: 




(\ - ^ sin 2 0d0 = AaE ^l-a 2 /b 2 



where E(k) = f* /2 \/ 1 - k 2 sin dO is the complete elliptic integral of the second 
kind. Thus the averaged equations for r± are simply: 

(12) f± = Af ± - f 3 ± c, 
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where c = 4a£( v / l - a 2 /b 2 ). 

The analysis of the MFI bifurcation is now straightforward. First consider the 
equation for f+. For all parameters, this has a hyperbolic attracting fixed point 
at the largest real solution of r 3 — Ar — c = 0. According to Theorem 4.1.1 in [8] 
the original equation has a hyperbolic attracting periodic orbit that passes close to 
the hyperbolic fixed point on a Poincare section. This periodic orbit 7 + encloses 
an MFI set. Further it is can be shown that any orbit beginning outside of 7 + will 
converge to the MFI set. 

Next consider the equation for f_. A saddle- node bifurcation occurs in the 
averaged equation when 

3c 2 / 3 

(13) Abif = 

at 



(14) r*=(f) : 



in which a stable and an unstable fixed point are formed. It follows from The- 
orem 4.3.1 in [8] that the Poincare map for original equation also undergoes a 
saddle-node bifurcation near the one for the averaged equation. Thus a stable peri- 
odic orbit 7~ for r_ is created. The stable orbits 7 + for r + and 7" for r_ enclose 
an annular MFI set. 

We summarize the above calculations. For A < Abif the disk region inside 7+ is 
a minimal forward invariant set. For A > Abif the MFI set is the annular region 
bounded by 7" and 7+. Thus the MFI set changes discontinuously at Abif and so 
a hard bifurcation of type B2 occurs. □ 

3. Random cycles 

The deterministic Hopf bifurcation involves the creation of a limit cycle. In this 
section we comment on the occurrence of attracting random cycles. 

Random cycles are defined in analogy with random fixed points [1]; for its def- 
inition we need to consider the noise realizations £ and the flow <j>^(x, £) for two 
sided time t G EL We henceforth consider the skew product flow 

with 

0*6 = 

for t, s G EL We formulate a result with a relaxed notion of attracting random 
cycle, as we allow for time reparameterizations. The notion implies though that a 
Poincare return map on a section {0 = 0} has an attracting random fixed point. 
Recall that a random fixed point is a map R :U — R 2 that is flow invariant, 

for P almost all £. A random cycle is defined as a continuous map S : U x S 1 — >• R 2 
that gives an embedding of a circle for P almost all £ G U and is flow invariant in 
the sense 

(s^s^ocLK^^s 1 ),?),*) 

Different regularities of the embeddings S 1 i->> ^(^S 1 ) may be considered. 
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The following result establishes the occurrence of attracting random cycles fol- 
lowing the hard bifurcation, for small noise amplitudes. We confine ourselves with 
a statement on continuous random cycles, thus avoiding for instance constructions 
of invariant cone fields to prove Lipschitz continuity. 

Theorem 3.1. For values of (A, e) with A > Abif and e small, the MFI E\ is 
diffeomorphic to an annulus and the flow $\ admits a random cycle S '.UxS* 1 —> R 2 
inside E\. 

The random cycle is attracting in the sense that there is a neighborhood U\ 
of the MFI E\, so that for all x G U\, the distance between (<£^(x,£),£) and 
\J te -^(S(0 t ^ S 1 ), t) goes to zero as t — >> oo. 

Proof. Recall from the proof of Theorem 2.1 the blow-up differential equations in 
polar coordinates: 

'=*-^/^> (£,')). 

(15) i = c -V3 + l/ A ( «W - si " 

i = 1. 

where u, v are noise realizations. 

Consider a reparameterization of time, r = g(f, u, v)t with 

9 ,m,,.,„=i +e ^(,(:),(-™; 

Note that g is close to 1, in particular positive, for small values of e. The reparam- 
eterization yields differential equations, for functions f, 0, t of r, 

-/ 1 ( \ — -3 . i/3^*/-n t / * f u \ ( cos(9 



(16) 



g(r,0,u,v) V 

-.-2/3 



g{r,6,u,v) 

Note that the time reparameterization does not preserve the flow of individual 
cycles S 1 ). However, for a fixed noise realization £, the graph |J tGR (5'(^ t ^, S 1 ), t) 
is an invariant manifold for both (15) and (16). As the time parameterization has 
changed only by a factor ^ f q u ^ that is bounded and bounded away from zero, to 
prove the existence of an attracting random cycle it suffices to do this for (16). 

The differential equations (16) define a skew product flow 

Define a graph transform T\ on embedded circles written as graphs U x [0, 2tt] \-> 
[R-,R+] for suitable < R- < R+. That is, T T X is determined through the property 

graph T\T% •) = *I(graph T{9~ T ^ ■), 0"^) 

Pick such that, in the limit e = 0, {R- < r < R+} is invariant. This 

is possible for A > Abif; note R- < r_ < r + < R+. Note that £ maps 
C°([0,2tt], [#_,#+]) into itself. We obtain S(£, •) as the limit 

(it) s(f,-)= iim rjr(c,-) 
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computed in the supnorm on C°([0,27r], 

Consider the translation s = f — R- . Then s satisfies the differential equation 

s' = (A - 3R 2 _)s - 3Rs 2 - s 3 + A#_ - i? 3 , + a, 
where we abbreviated a = /a( ^ ) , I C ° S ? I \ • We consider the differential 



v I \ sm( 

equation for s on [0,i?+ — R-], which is invariant by the choices of Note 
also that the coefficient A — 3R 2 _ is negative. Suppose s\ > S2 are two solutions and 
consider the difference u = s± — #2- Then u satisfies 

v! = (A - 3R 2 _ - 6R-s 2 - 3s 2 2 )u - (3R- + 3s 2 )u 2 - u 3 . 

All coefficients here are negative, implying that u(r) —> as r — >> oo. This compu- 
tation demonstrates the convergence in (17). Likewise for small values of e. □ 

Write f H> H\(r, ^) for the first return map on {# = 0} of 3>t(f, 0, <f). As a 
corollary of the above theorem we obtain an attracting random fixed point for IIa- 

Corollary 1. For A and e as in Theorem 3.1, the first return map Ii\ admits a 
random fixed point R{C) inside E\ D {0 = 0}. 

The random fixed point is attracting in the sense that there is a neighborhood U\ 
of E x , so that for all x eU x n{6 = 0} ; \U k (x^) - n k (R(£))\ as k oo. 

Proof Just note that IIa = ^ r ^ 7r£2/3 , defined in the proof of Theorem 3.1. □ 

A better notion of random attracting cycle would be without time reparameter- 
izations; 

$i(5(€,S 1 ),0 = 5(^ ) S 1 ). 

This would allow a discussion of the dynamics on the random cycle, such as the 
possibility to find an attracting random fixed point on it, compare [2, 3]. 

4. Simulations of bounded noise, radially symmetric Hopf. 

4.1. Radially symmetric Hopf with bounded noise. In this section we repro- 
duce invariant measures for the radially symmetric system 

x = Xx — y — x(x 2 + y 2 ) + eu, 

y = x + Xy — y(x 2 + y 2 ) + ev, 

where A G R is a parameter and u and v are noise terms satisfying: u 2 + v 2 < 1, 
representing radially symmetric noise. In the simulations below we generate the 
noise by stochastic differential equations with reflective boundary conditions. 

For e = the differential equations (18) undergo a supercritical Hopf-Andronov 
bifurcation at A = 0. For A < the origin is a stable global attr actor. For A > 
the origin is unstable, but there is a circular globally attracting periodic orbit at 
r = vX Changing to polar coordinates we find: 

r = Xr — r 3 + sa, 

(19) . £ 

= 1 + -(3, 
r 

where 

a = (u cos + v sin 6) = (u, v) • (cos sin 0) 

and 

j3 = {—u sin + v cos 0) = (v, —u) • (cos ^, sin 6) . 
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Figure 1 indicates boundaries of radial components of the MFI set for (19). Follow- 
ing the proof of Theorem 2.1 the boundaries of the MFI set are given by zeros of 
the upper and lower radial differential equations 

(20) f± = \r± - r4 ± e. 

From (13) and (14), 

Write p + for the positive zero of the upper differential equation for r+. For 
A > Abif, let p- be the largest positive zero of the lower differential equation for 
r_. For A < Abif the MFI set is a disk of radius p + . For A > Abif the MFI set is an 
annulus bounded by circles with radii See Figure 1. 




Figure 1. Boundaries of the radial components of the minimal 
invariant set in the random symmetric Hopf bifurcation (18). Here 
£ = 0.05. The hard bifurcation takes place at Abif ~ 0.2565 and 
the inner boundary has initial radius r* « 0.2924. The MFI set 
changes from a disk with radius p + for A < Abif to an annulus with 
radii p_ < p + for A > Abif. The graph labeled r s indicates the 
boundary of the set of fixed points which are stable for A < and 
unstable for A > 0. 

In addition to the MFI set, another important dynamical feature of these equa- 
tions is the set of equilibrium points. For A < these equilibrium points are stable, 
but for A > they are unstable. These are the points that are fixed points of (18) 
when (u, v) is a fixed value. One can find that this set is a disk centered at the 
origin with radius r s which is a solution of: 

(l-2A)r 6 + (l + A 2 )r 2 -£ 2 =0. 

The radius r s as function of A is plotted in Figure 1. 
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4.2. Simulations. In this section we simulate the bounded noise bifurcation for 
the sake of demonstration. We note that associated with MFI sets are invariant 
densities whose supports are the MFI sets [9] . In the symmetric Hopf bifurcation 
as in the previous subsection, we approximate these invariant densities for different 
values of A as A moves through the random bifurcation. 

The generation of bounded noise is somewhat arbitrary without specific knowl- 
edge about the noise involved in a particular setting. In this work we generate the 
noise terms u and v via the stochastic system: 

du = (7 dW\ , 

( 21 ) dv = adW 2l 

where dW\ and dW 2 are independent (of each other) normalized white noise pro- 
cesses and (21) are interpreted in the usual way as Ito integral equations. In order 
to assure boundedness and radial symmetry, we restrict (u, v) to the unit disk and 
impose reflective boundary conditions. Other methods of generating and bounding 
the noise did not produce significantly difference in the results. 

The parameter a can be interpreted as the rapidity of the noise. If a is small, then 
u and v change slowly and as a increases, they change more quickly. It turns out 
in our simulations that the value of a has a strong influence on the characteristics 
of the invariant density of (18). 

In the first set of simulations Figure 2, a = 1.0 as an example of fast noise and 
in Figure 3, a = 0.00001 to show slow noise. Note that for slow noise and A > 
there is a separation of time scales. Specifically, for each A and each set of values of 
(u, v) the deterministic system defined by holding (u, v) fixed has an exponentially 
attracting periodic orbit. We expect then that the invariant density is concentrated 
on the set which is the union of all of these deterministic limit cycles. For fast noise, 
this will not be the case and as is observed in Figure 2, the approximated invariant 
densities tend to be smoother than for slow noise. 

In both sets of simulations we use e — .1. This leads to the following values of 
the bifurcation parameter and the radius of the inner boundary of the MFI set at 
the bifurcation: 

Abif = ~ 0.407163, r* = (-) 3 ^ 0.368403. 

We began the simulations by selecting random initial values for (u, v) and (x, y). 
As we did not need an accurate solution to generate the noise we approximated 
solutions to the SDE in (21) using a first order Taylor method, namely Euler's 
method. We then added this noise into the system in (18), which we solved using 
a second order Adam's-Bashforth method. 

After selecting the initial noise term and a starting point for the RDE, we ran 
the system for 1000 iterations to allow the solution to the RDE to move into the 
MFI sets. We then ran the system for 5 x 10 5 iterations and recorded the values of 
every fifth (x, y) coordinate. We repeated this process for 10 2 starting points and 
recorded a total of 10 7 points, which we used to generate a 2-dimensional histogram. 
Bright regions indicate that large numbers of samples were observed in this region 
indicating a higher value of the invariant density. Darker regions indicate lower 
density. 

The deterministic bifurcation takes place at A = and the random bifurcation 
occurs at Abif ~ 0.4072. In each set of figures, the densities are plotted for multiples 
of Abif; namely for A = .01 Abif, .05 Abif, .1 Abif, .5 Abif, Abif, and, 1.1 Abif. The outer 
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Figure 2. Invariant densities for fast noise: a = 1 and e = 0.1. 




Figure 3. Invariant densities for slow noise: a = .00001 and e — 0.1. 

circle is the outer boundary of the MFI set, the circle inside the MFI is the stable 
periodic orbit of the deterministic system and the inner circle that appears at Abif 
is the inner boundary of the MFI set. It has initial radius r* = 0.3684. 

The invariant densities are positive for points inside the MFI set, but can be 
expected to go to zero very rapidly as the boundary is approached (compare [14]). 
Also note that the densities become undetectably small in the center long before 
the random bifurcation. 
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A color movie of the simulations can be found at: 
http://www.youtube.com/watch?v=4tVtWGdVMi8. 
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